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These are personal recollections of how Gerry Brown's work and thinking 
■ influenced the development of the chiral effective Lagrangian based theory 

CNJ ! of nuclear forces. 

> : 

■ 

in ; 1 The early days at Stony Brook 

I first met Gerry on the occasion of the Erice school in 1981, where he and Mannque 
' Rho were giving quite interesting lectures on the chiral bag model and related issues. At 

that point, I was close to finishing my diploma thesis under the supervision of Manfred 
Gari at Bochum, working on two-meson exchange currents in the framework of unitary 

1 ^ i ' transformations. More specifically, the goal was to find out the relative contribution of 

the isoscalar p7r 2 7 current compared to the well-established leading order isoscalar pir^y 
current. This calculation turned out to be - to use Howard Georgi's words - an excellent 

£T) ' exercise in self-torture, as one had to evaluate about 1000 time-ordered diagrams by 

hand. Also, the result that these corrections were already sizeable at momentum trans- 
fers below q 2 ~ 1 GeV 2 did not quite fit into the standard lore which was so eloquently 
summarized by Gerry and Mannque in their chiral filter hypothesis pQ. This hypothesis 

t— i ■ that grew out of explicit meson-exchange current calculations including heavier hadrons 

like the A(1232) or the vector mesons p, cj was clearly emphasizing the role of the pion. 



Whenever present, the one-pion exchange current was supposed to dominate the perti- 
nent response of light nuclei to electromagnetic probes, even up to momentum transfers 
as large as 1 GeV 2 , way beyond the soft-pion limit. In fact, it was clear to Gerry and 
Mannque - and others - that the understanding of the nucleon structure from QCD 
and the nuclear interactions were two intimately connected issues, firmly rooted in the 
chiral symmetry of QCD. It is worth quoting from Ref. [Ij: "... that there must exist an 
intimate connection between what makes up a hadron and what induces a rich variety 
of interactions between hadrons. it will be rather unlikely one will gain a full un- 
derstanding of one without the other." Only then too much emphasis was put onto the 
quark structure of the nucleon - as we understand now, low-energy QCD is well approx- 
imated by an hadronic effective Lagrangian. When doing my calculations at Bochum, I 
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did not know that, but learned about it quickly when I arrived at Stony Brook in April 
of 1982. The Erice lectures on the little bag model had aroused my interest and, hon- 
estly, I was fed up calculating meson-exchange currents. Therefore, I contacted Gerry 
about becoming a graduate student at Stony Brook. It did not take him long to accept 
me - as it turned out, I spent a short but very productive time in Gerry's group until 
my graduation in December 1984. But let me briefly come back to the meson-exchange 
current (MEC) story. In fact, only after some strong support from Dan-Olof Riska, one 
of the MEC pioneers [2J, my diploma work was finally published [3] but largely ignored 
by the community. The crux of the matter was that in those days one did not have a 
power counting that allowed one to systematically address the corrections to the leading 
order one-pion contribution either in the nuclear forces or in the corresponding exchange 
currents, though a tremendous amount of phenomenology had been built up over the 
years. In particular, the Stony Brook and Paris approach to the two-nucleon force that 
used dispersion relations to connect the irN scattering amplitudes to the 2-7r-exchange 
in the NN force supplemented by some heavy mass Yukawa-type functions with param- 
eters determined from a fit to the data where already very close to what would become 
later the "modern theory of nuclear forces". In any case, I had to learn the basis of the 
fundamentals of the Stony Brook potential by working through the book of Gerry and 
Andy Jackson [4J. This turned out to be a very painful exercise as the notation would 
change from chapter to chapter, but made me familiar with dispersion theory and the 
many intricacies of nuclear interactions. In fact, one of the early projects Gerry assigned 
to me was the coupling of p-mesons to the little bag, which I worked out within a few 
weeks and which was supposed to be the starting point for a more microscopic investi- 
gation of vector meson exchanges in the nuclear force, thus further contributing to the 
sought-after unified picture of nucleons and the nuclear interactions. When I presented 
to Gerry my notes in the form of a handwritten draft - ETgXwas not yet available - 
he just told me that I was too fast for him and proposed to me that I work alone. So 
I teamed up with Andreas Wirzba - a fresh new graduate student from Miinster - and 
Ismail Zahed, then a post-doc in Gerry's group to work on Casimir effects in chiral bag 
models and on many aspects of Skyrmions, which was a hot topic at that time. The 
Skyrme model promised to lead to the long wanted unification of nucleon and nuclear 
physics, and many impressive calculations revealed further deep connections, but again 
the lack of a power counting was the stumbling stone. 



2 A first step: Studying pion-nucleon scattering in chiral 
perturbation theory 

During my stay as a Heisenberg fellow at Bern, I was lucky to learn chiral perturbation 
theory (CHPT) from Jiirg Gasser and Heiri Leutwyler, who had transformed Weinberg's 
ideas into an effective machinery to analyze low-energy QCD. With Veronique Bernard 
and Norbert Kaiser we systematically analyzed electroweak reactions on the nucleon 
and in particular neutral pion photoproduction off the proton, — > ir°p, triggered 
by measurements at Saclay, Mainz and Saskatoon to test the then believed low-energy 
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Figure 1: The LECs Cj (circles) in pion-nucleon scattering (left), the two-nucleon (NN) 
interaction (center) and the three-nucleon (NNN) interaction (right). 

theorems based on soft-pion algebra (plus one extra assumption) . The explanation of the 
data by a new one-loop effect helped to establish the method as a precision tool in the 
single nucleon sector. In the the late nineties and early years of the new millennium, with 
my students Sven Steininger and Nadia Fettes [5] we investigated pion-nucleon scattering 
in the framework of heavy baryon chiral perturbation theory - this later turned out to be 
an important ingredient in the construction of the few-nucleon forces based on effective 
field theory (EFT). The underlying chiral Lagrangian of pions and nucleons coupled to 
external currents can be written as 

C cS = + + £0) + £(4) + . . . , (1) 

where the superscript denotes the chiral dimension (the power in small external momenta 
and/or pion masses). Tree level computations are done with + C^ 2 \ a complete 
one-loop calculation involves all loop graphs with insertions from jC^ 1 ' and at most one 
insertion from and so on. In its full glory, £ fj can be found in [7]. Of particular 
importance in these studies are the finite dimension two low-energy constants (LECs) c,, 
see Fig. [H that also feature prominently in the description of the two- and three-nucleon 
forces - Gerry was one of the first who had understood this long before the advent 
of EFT. An even better determination of some of the q was possible by analytically 
continuing the amplitudes into the interior of the Mandelstam triangle, see Fig. [2] - the 
ability to use dispersion relations as I learned it in Stony Brook turned out to be very 
valuable. All this work can be summarized by the following values for the Cj in units of 
GeV- 1 : 

a = — 0.9±g;| , c 2 = 3.3 ± 0.2 , c 3 = -4.7±i;g , c 4 = 3.5±g;i . (2) 

The values for C2, C3 and C4 are larger than the expected (natural) values: q ~ 5a/A x ~ 
1.1 . . . 1.5, with gA = 1.267 the nucleon axial-vector coupling and A x ~ 1 GeV the 
scale of chiral symmetry breaking. This was understood early in terms of resonance 
saturation - the fact that the heavier states integrated out from the EFT leave their 
imprint in the values of the LECs. This method was pioneered in the meson sector 
by Gerhard Ecker and collaborators and also by John Donoghue and collaborators [5]. 
In 1997, with Veronique Bernard and Norbert Kaiser, we had already determined the 
Cj by a comparison to scattering lengths and subthreshold parameters and extended 
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u = (m+ M) 2 \ 



Figure 2: Mandelstam plane. The Mandelstam triangle is the inside of the thick lines. 

the resonance saturation scheme to the one-baryon sector, see Fig. [3J Not surprisingly, 
the LECs C2 and C3 are largely dominated by the close-by A(1232) resonance, whereas 
C4 also receives an important contribution from the p-meson |9j. This leads to the 
following values of 02,3,4 (in brackets the ranges obtained from resonance saturation) C2 = 
3.9 [2 . . . 4], C3 = —5.3 [—4.5 ... — 5.3], C4 = 3.7 [3.1 . . . 3.7], whereas c\ is given by scalar 
meson exchange (tttt correlations). These numbers are also consistent with the recent 
analysis of pion-nucleon scattering and subthreshold parameters derived from pionic 
hydrogen and deuterium, c\ = — 1.2 . . . — 0.9, ci = 2.6... 4.0, C3 = —6.1... — 4.4 [TU] . 
Let me also stress that the uncertainties in the Cj are strongly correlated, as best seen 
from the combination — 2c\ + C2 + C3 
pion-nucleon scattering length. 



g\/(8m) that determines the small isoscalar 
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Figure 3: Resonance saturation for the dimension two (circles) and three (squares) LECs. 



3 Flashback: The Texas revolution 

Already in the early nineties a decisive step towards a systematic theory of the nuclear 
forces had been done. This was initiated by a series of papers by Steven Weinberg, who 
had already pioneered the use of effective chiral Lagrangians decades earlier under the 
nowadays somewhat strange-reading title "Phenomenological Lagrangians" [TI]. In that 
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marvelous paper the idea of the power counting was developed in detail and further a very 
illuminating application of the renormalization group to work out the so-called chiral 
logarithms was presented. Years later, while teaching the effective field theory of pions 
and nucleons, Weinberg suddenly realized that there are additional terms consistent 



with the power counting, terms quadrilinear in the nucleon fielda^H He immediately 
realized that these must be related to the much discussed mechanism for the short-range 
repulsion of the nuclear forces. Here the fundamental principle of effective field theory 
works at it best - in the low-energy sector of a given theory one is not able to resolve 
the physics at large momentum scales but simply parameterizes it in terms of operators 
made of the fields one has at one's disposal - it just does not matter whether the nucleon- 
nucleon repulsion is due to vector-meson exchange, quark rearrangement energies or 
whatever model has been invented before. Weinberg then worked out the consequences 
of the spontaneous chiral symmetry breaking for the forces between two, three and four 
nucleons using the method of phenomenological Lagrangians |12] , nowadays called chiral 
perturbation theory or chiral effective field theory. Remarkably, the first of these two 
papers contains only three references, one to the review Gerry had written in 1985 with 
Sven-Olaf Backman and Jouni Niskanen [13| and also, he explicitly thanks Gerry for 
"enlightening conversations on nuclear forces". The important step taken by Weinberg 
was to realize that in the nuclear force problem one can not apply the power counting 
directly to the S-matrix, but rather to the effective potential - these are all diagrams 
without ./V-nucleon intermediate states. Such diagrams lead to pinch singularities in the 
infinite nucleon mass limit (the so-called static limit), so that e.g. the nucleon box graph 
is enhanced as m/Q 2 , with m the nucleon mass and Q <C A x a small momentum. The 
beautiful power counting formula for the graphs contributing with the uth power of Q 
or a pion mass to the effective potential reads (considering only connected pieces): 



2 - N - 2L + V i 



di H 

1 2 



(3) 



Here, N is the number of in-coming and out-going nucleons, L the number of pion loops, 
Vi counts the vertices of type i with dj derivatives and/or pion mass insertions and n% is 
the number of nucleons participating in this kind of vertex. Because of chiral symmetry, 
the term in the square brackets is larger than or equal to zero and thus the leading terms 
contributing e.g. to the NN potential can easily be identified. These are the time-honored 
one-pion exchange and two four- nucleon contact interactions without derivatives. These 
contact interactions were indeed the missing part which were earlier modeled by heavy 
meson-exchanges or by some kind of fit function - but none of these approaches was 
controlled or systematic. Here again the awesome power of power counting becomes 
crystal-clear - it just took one person not biased by nuclear folklore to take this final step. 
The so-constructed effective potential is then iterated in the Schrodinger or Lippman- 
Schwinger equation, generating the shallow nuclear bound states as well as scattering 
states. The resulting contributions to the 2N, the 3N and the 4N forces are depicted 
in Fig. [H I will come back to a detailed discussion of the various entries below, at this 



#1 I am grateful to Steve Weinberg for providing me with this recollection. 
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Figure 4: Contributions to the effective potential of the 2N, 3N and 4N forces based 
on Weinberg's power counting. Here, LO denotes leading order, NLO next-to- 
leading order and so on. Dimension one, two and three pion-nuclcon interactions 
are denoted by small circles, big circles and filled boxes, respectively. In the 4N 
contact terms, the filled and open box denote two- and four-derivative operators, 
respectively. 

point it is important to observe that - consistent with phenomenological observations 
- three-nucleon forces appear only two orders after the dominant NN forces and four- 
nucleon forces are even further suppressed, appearing only at N 3 LO. Let me stress that 
the beautiful analysis of few-nucleon forces in chiral EFT by Bira van Kolck |14| was 
an important ingredient to set up the power counting as displayed in Fig. [J] and that 
this extended power counting has been the major development on the way to the chiral 
theory of nuclear forces. 

The first large-scale numerical evaluation of this chiral potential was done by Bira van 
Kolck and collaborators a few years after Weinberg's seminal papers, see Refs. |15j . They 
worked up-to-and-including N 2 LO in a theory of pions, nucleons and also the A isobar, 
fitting the LECs of the contact interactions to low-energy NN phases separately for 
isospin 1 = and 1 = 1 and properties of the deuteron. This work was certainly ground- 
breaking but it was also not immediately accepted in the nuclear physics community - 
for the simple reason that the phenomenological approaches to the two-nucleon problem 
had achieved a much higher precision - at the expense of many more parameters and lack 
of systematicity. Still, the first important step was done. Before picking up this theme, 
let me stress that Weinberg extended his approach to pion reactions on light nuclei |16| . 
using at that time deuteron wave functions from phenomenological potentials. This so- 
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called "hybrid appraoch" laid the ground for systematic studies of pion scattering off 
light nuclei, pion photo- and electroproduction off light nuclei and of pion production in 
NN collisions. 

4 The chiral theory of nuclear forces as a precision tool 
4.1 The nucleon-nucleon interaction 

When I was still at Bonn filling in for the vacant chair of Max Huber, then rector of Bonn 
University, I met a few times with Karl Holinde of the IKP Jiilich. Karl had become 
famous for his work on the Bonn potential and was hired by Josef Speth at Jiilich in 
1984, which gave Karl a fantastic home base to continue his work. We talked about the 
chiral potential of the Texas group, and he was clearly interested but also aware that 
the precision was not yet competitive. So we decided to work together to try to improve 
the Texas potential. Unfortunately, before really getting started, this collaboration was 
terminated by Karl's untimely death in December of 1996. In fact, Gerry came to 
Jiilich the next year and gave a wonderful talk about the roots, the physics and the 
achievements of the Bonn potential, focusing on Karl's contribution. Fortunately, a bit 
later Walter Glockle from Bochum contacted me - he had acquired a brilliant young 
student named Evgeny Epelbaum, who had done some calculations of the chiral NN 
potential in his diploma thesis and was supposed to continue working on this topic 
for his doctoral thesis. Walter was (and still is) a world-leading expert in few-nucleon 
calculations but did not feel at ease with effective chiral Lagrangians. So we decided 
that we would team up to supervise Epelbaum, which turned out to be a true pleasure 
as we did not have to do much. Being aware of the short-comings of energy-dependent 
potentials in few-body calculations, we decided to recalculate the Texas potential using 
the method of unitary transformations - which I had not used for about 15 years. Based 
on Epelbaum's diploma thesis, we published a paper on the power counting adapted 
to this method in [T7j. Before working out the phenomenological consequences of this, 
we decided to work out some simpler examples to get better acquainted with the many 
subtleties of the nuclear interaction and how they reflect on the corresponding effective 
theory. Interestingly, in two papers we developed a method to integrate out the high- 
momentum components from an effective theory based on a simple local potential of the 
Malfiiet-Tjon type [IB]. This preceeded the work of Achim Schwenk, Tom Kuo, Bengt 
Friman, Gerry and others at Stony Brook on the so-called Mow-k potentials based on 
elegant renormalization group techniques, see e.g. the review |19] and references therein. 
Such parallel developments are a fine witness of Gerry's legacy - he has taught many 
generations of students a deep understanding of physics that will necessarily lead to 
progress in the field of nuclear physics. But back to the nuclear force problem. We 
published our N 2 LO results in [2D], which were markedly improved as compared to 
the Texas potential. The simple reason was that we did not make global fits but rather 
projected the chiral potential onto partial waves and determined the four-nucleons LECs 
by fitting to the S- and P-waves and the 3 Si- 3 -Di mixing parameter in the neutron-proton 
system. An important paper that had been published earlier was the work by the Munich 
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Figure 5: Left: np S-waves at NLO (grid), N 2 LO (light shaded) and N 3 LO (dark shaded). 

Right: Differential cross section and vector analyzing power for np scattering at 
E U h = 50 MeV compared to the N 2 LO (light shaded) and N 3 LO (dark shaded) 
predictions. 



group, where the role of the two-pion exchange to the peripheral partial waves was 
scrutinized [2T]. They demonstrated that these peripheral waves are well described by 
one- and two-pion exchanges in the Born approximation, thus only the low partial waves 
were to be used in fitting the multi-nucleon LECs. However, the large values of the Cj, cf. 
Eq. (J2J) > lead to a very strong isoscalar central potential at short distances that generates 
unphysical deep bound states for the larger values of the cut-off scale (A ~ 1 GeV) in the 
LS-equation. Again, earlier work done by Chemtob, Durso and Riska - that was taught 
at Stony Brook - came to our rescue, we simply adopted the so-called spectral function 
regularization method to deal with the unwanted short-distance behavior of the two-pion 
exchange [22J. So the ground was prepared for going to yet higher orders. While Evgeny 
had moved to Jefferson Lab as the first Nathan- Isgur-Distinguished-Fellow, we continued 
our collaboration and published the N 3 LO results in [23]. Earlier on, Ruprecht Machleidt 
and David Entem had published N 3 LO results using dimensional regularization in the 
two-pion exchange [24] based on a series of papers from Norbert Kaiser - this story is 
told in the contribution of Machleidt to this volume. Some characteristic results of our 
N 3 LO calculation are displayed in Fig. [5] - these are not more accurate than the ones 
based on phenomenological potentials, however, they provide a measure of the theoretical 
uncertainty and, more importantly, they are based on the chiral symmetry of QCD. To 
close the circle, one can indeed show that the four-nucleon operators to a large extent 
are saturated by resonance excitations, as explicit calculations mapping boson-exchange 
models onto the chiral expansion of the effective potential showed, for details see |25| . 
In a way, EFT provides a reason why the meson-exchange models of the nuclear force 
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Figure 6: Topologies of the leading contributions to the chiral 3NF. From left to right: 
Two-pion exchange, one-pion-exchange and 6N contact interaction. 

have been so successful. 

4.2 Three-nucleon interactions 

One of the biggest advantages of the chiral Lagrangian approach to the nuclear force 
problem - as stressed early in Weinberg's papers - is the consistent derivation of the 
3- and 4-nucleon forces and also the meson-exchange currents (that is, the response to 
electroweak probes). It is worth mentioning that decades earlier Gerry was aware of the 
role of chiral symmetry in the description of three-body forces in nuclei - see his nice 
paper with Saul Barshay in 1972 [25] • But let me come back to the modern approach. As 
already noted, the three-nucleon force (3NF) only appears two orders after the leading 
NN interaction. At this order, there are only three topologies contributing, see Fig. O 
The two-pion exchange topology is given again in terms of the Cj, as discussed in detail 
in [27] . The so-called D-term, which is related to the one-pion exchange between a 4N 
contact term and a further nucleon, has gained some prominence over the last years, 
as many authors have tried to pin it down based on a cornucopia of reactions, such 
as Nd -> Nd, NN -> NNtt, NN ->■ d£v e or dir jNN. This demonstrates nicely 
the power of EFT - very different processes are related through the same LECs thus 
providing many different tests of chiral symmetry (as it is also the case with the LECs 
Cj, see Fig. [T]). The LEC E related to the 6N contact interaction can only be fixed in 
systems with at least three nucleons, say from the triton binding energy. Remarkably, 
although there are many diagrams at N 3 LO, there are no new LECs to be determined. 
The long-range terms of this force can be found in [25] and the shorter ranged ones will 
be published soon (see also [2J5] ). Many applications of these forces and the testing of 
their structure can be found in the reviews [3(J|. I3T] . 

As a very nice example that these forces make their way into nuclear structure cal- 
culations, I show here the results from Petr Navratil and collaborators [32], who per- 
formed large-basis no-core shell model (NCSM) calculations including the leading chi- 
ral 3N forces and demonstrated their necessity to describe the spectra of nuclei with 
A = 10, ... 13, see the entry NN+NNN in Fig. [JJ Note, however, that in this calculation 
the 2N force was employed at N 3 LO while the 3NF was only included at N 2 LO. It will 
be interesting to see how the inclusion of the sub-leading 3NF terms will modify these 
results. 

So far, I have only discussed the EFT with pions and nucleons. Clearly, the two-pion 
exchange diagram in Fig. \5\ is a controlled approximation to the time-honored Fujita- 
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Figure 7 : Spectra of light nuclei using chiral NN and chiral NN+NNN forces in a NCSM 
calculation compared to the data. Figure adopted from [32] and courtesy of 
P. Navratil. 



Miyazawa force, whose 50 th anniversary was celebrated in Tokyo in 2007. For a light 
reading of how the Fujita-Miyazawa force is related to the EFT description of the 3NF, 
I refer to [33] and references therein. 

4.3 The nuclear matter problem 

Another important question in nuclear physics is to understand the saturation properties 
of nuclear matter - an idealized infinite system of nucleons in which all Coulomb effects 
are switched off. From the properties of heavy nuclei using some sophisticated mass for- 
mula, one can extrapolate to nuclear matter - and determine its saturation properties. 
The binding energy per nucleon E/A is approximately —16 MeV at a Fermi momen- 
tum of about 1.3 fm . So what does chiral EFT have to say? A first important step 
was taken by Norbert Kaiser, Wolfram Weise and collaborators at Miinchen, who calcu- 
lated the contribution of pion exchange(s) to the energy density of nuclear matter and 
showed that the energy density of isospin symmetric nuclear matter can be extremely 
well approximated by the simple form 

q/>2 l3 L4 

10m wr m J 

They then calculated the coefficients a and j3 from chiral dynamics. With some fine- 
tuning of the regulator, one finds an astonishingly good description of the energy density 
of nuclear matter [34] . This was later improved by including e.g. higher orders (sensitive 
again to the LECs q) and isobar degrees of freedom [35] . Interestingly, all this was based 
on a loop expansion with no explicit power counting - a power counting that indeed 
explained the success of these calculations was only set up much later, see [3S]- Space 
forbids a discussion of this power counting and the resulting physics in detail - I only 
would like to mention that it could be shown that for many reactions the contributions 



10 




k F [fm '] k p [fin" 1 ] 



Figure 8: Left panel: Binding energy per particle of nuclear matter as a function of the 
Fermi momentum kp for the N 2 LO potential using a Lippman-Schwinger cut off 
A = 550 MeV and a spectral function cut off A = 600 MeV (solid line). The 
binding energies obtained in the absence of three-nucleon interactions are shown 
by the dashed line. The black square gives the empirical nuclear matter prop- 
erties. Right panel: Saturation points of nuclear matter. Downward triangle, 
upward triangle and rectangle: various phenomcnological approaches. Circles: 
chiral EFT at N 2 LO. 



from the multi-nucleon interactions cancel to leading order which is at the heart of the 
success of the Munich group calculations. 

Here, I briefly pursue a more "conventional" approach to the nuclear matter problem 
based on more standard nuclear many-body theory. For decades, the method of choice 
to analyze nuclear matter has been the G-matrix of Brueckner and others - leading 
to the problem that for phase-shift equivalent NN potentials the resulting saturation 
properties lie on the Coester line that does not pass the empirical values. Recently, 
Siegfried Krewald and collaborators have recalculated the properties of nuclear matter 
based on the chiral NN and 3N forces (see the v2 version of Ref. |38j ) utilizing the R- 
matrix as proposed by Baker in 1971 [37J. As can be seen from the left panel of Fig. [51 
the NN interaction alone only binds for Fermi momenta larger than 1.7 fm -1 . However, 
including the 3NF with natural values for D and E, one finds binding with the proper 
strength at the proper density. In the right panel, the sensitivity to the two cut-offs in the 
NN force as compared to phenomenological determinations is shown, it is of comparable 
size. The dependence on the 3-body LECs D and E is also weak, see Fig. 3 in Ref. |38| . 
It will be interesting to push these calculations to N 3 LO and to include the A(1232) in 
the EFT. 

4.4 Back to square one - meson-exchange currents 

As told in section 1, 1 started my career calculating meson-exchange currents based on the 
method of unitary transformations. It is quite a nice turn of events that in chiral nuclear 
EFT one is now able to do this in a controlled and systematic manner, thanks to the 
power counting. In our approach, Stefan Rolling, Hermann Krebs and Evgeny Epelbaum 
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Figure 9: Predictions for the n-d and n- 3 He capture reactions. Chiral EFT: the curve 
labeled N3LO(S-L) includes pion loop corrections contact terms, while the curve 
labeled N3LO(LECs) includes in addition contributions from (higher order) pion 
exchanges and non-minimal contact currents. Pre-EFT approaches: SNPA and 
SNPA* denote up-to-date conventional calculations. Figure courtesy of Rocco 
Schiavilla. 
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Accessible by 
Lattice QCD 




Figure 10: Left panel: Euclidean space-time lattice with point-like nucleons on the lattice 
sites and a typical lattice spacing of a = 2 fm. Right panel: Nuclear phase dia- 
gram as accessible by lattice QCD and by nuclear lattice EFT. Figure courtesy 
of Dean Lee. 



have taken up the task - the two-pion exchange electromagnetic current based on this 
method is published in and together with the Cracow group led by Henryk Witala 
and Jacek Golak a thorough investigation of electro-nuclear processes is underway, for 
a first application to deuteron photodisintegration see [3D]. However, the group around 
Rocco Schiavilla, using old-fashioned time-ordered perturbation theory, has taken the 
lead and performed calculations of the MECs of one- and two-pion range and applied this 
to magnetic moments of light nuclei [UJ. Most recently, they have calculated thermal 
neutron capture on d and 3 He |42j . The resulting radiative capture cross sections a^ d 
and o r j i 3 He and the photon circular polarization parameter R c , resulting from the capture 
of polarized neutrons on the deuteron, are shown in Fig. in comparison to two state-of- 
the-art conventional approaches. As the authors note, these processes are not the best 
or simplest to illustrate the convergence pattern of chiral nuclear EFT. So much more 
work is needed here and will be done. 

5 Flash forward: Combining nuclear EFT with lattice 
simulations 

Having now constructed nuclear forces between two, three and four nucleons, one has a 
new tool to address the nuclear many-body problem. One venue is to combine standard 
many-body techniques with these forces, as exemplified by the NCSM calculation shown 



13 



A 

UJ 

V 



-60 
-70 
-80 
-90 
-100 
-110 
-120 



LO 
+ ANLO 
+ AIB + AEM 
+ ANNLO 
+ 4 N Contact 

i -• - - 



20 
15 
10 
5 

-5 
-10 



1 1 1 

ANLO 
AIB + AEM 
ANNLO 
4/V Contact 


1 1 _ 

1 X ! 

!■ 

i — e — i 


""' D - a -ir4— 








1 1 1 





0.02 0.04 0.06 0.08 0.1 
f(MeV~') 



0.02 0.04 0.06 0.08 0.1 
f(MeV _1 ) 



Figure 11: Ground-state energy of 12 C as a function of Euclidean time. For details, 
see [35] . 



in Sec. 14.21 or the discussion of nuclear matter in Sec. 14.31 (see also the contribution by 
Josef Speth, Siegfried Krewald and Frank Grummer to this Festschrift). Here, I briefly 
want to outline a very different and novel approach that combines the chiral EFT for the 
forces with the method of Monte Carlo simulations, that are so successfully utilized in the 
lattice approach to QCD. The basic idea is to formulate the chiral EFT on a Euclidean 
space-time lattice as depicted in the left panel of Fig. [10]- here, the lattice spacing serves 
as an UV regulator and has to be chosen such that the nucleons - which are treated as 
non-relativistic point-like particles on the lattice sites - do not overlap. A lattice spacing 
of a = 2 fm entails an UV cut-off A = ir/a ~ 300 MeV. The pion and nucleon propagators 
as well as the one- and two-pion exchanges and the multi-nucleon contact interactions 
are written in terms of lattice variables, making use of Hubbard-Stratonovich fields to 
bring the 4- and 6- nucleon interactions into a quadratic form (for more details, see |43j 
and references therein as well as the nice review by Dean Lee |44] ) . Given this Euclidean 
formulation of the action, the generating functional of the theory can then be evaluated 
with stochastic methods, such as the hybrid Monte Carlo approach. It is important to 
stress that the approximate Wigner SU(4) symmetry of the nuclear interactions strongly 
suppresses the sign oscillations and thus makes such nuclear lattice simulations easier to 
handle than lattice QCD calculations. Within this scheme, a systematic study of nuclei 
up to A ~ 40 will be possible using petascale computing as long as the total spin and 
isospin of the nucleus under consideration is zero. Also, the nuclear phase diagram can 
be studied for a wide range of temperatures and densities, as indicated in the right panel 
of Fig. fTUl The state-of-the-art of these simulations is the calculation of the ground-state 
energy of 6 Li and 12 C [JS|. In these works, for the first time the contributions from the 
Coulomb interaction between protons and strong isospin-breaking effects were included. 
A parameter- free prediction for the energy difference between the triton and 3 He could 
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be given, 



£( 3 H) - £( 3 He) = 0.78(5) MeV , (5) 

in good agreement with the empirical value of 0.76 MeV. The ground-state energies 
of 6 Li and 12 C are calculated as -32.9(9) MeV and -99(2) MeV, respectively, not far 
from the empirical values of —32.0 and —92.2 MeV, cf. also Fig. [TTJ This accuracy is 
comparable to other so-called ab initio calculations (like the NCSM or Greens function 
Monte Carlo), that often are based on a less consistent formulation of the underlying 
nuclear forces. 

Taking 12 C as the benchmark, it is interesting to estimate the required computing time 
and storage needed for the calculation of larger nuclei^ 2 !. With an improved algorithm, 
the required CPU time for the 12 C simulation on a BlueGene/P architecture and the 
necessary storage to make the configurations (about 3500) available is 

X? 2 p c y = 5 x 10~ 5 PFlop - yr , xg£ age = 0.14 TB . (6) 

For a nucleus with A nucleons and with total spin S and isospin I the required CPU 
time and data storage place can be estimated as 

X CPU « XW^ x ^) 3 2 exp[0.10(^ - 12) + 3(5 mod 2) + 4/] , 

^storage _ xH°™ sc x f^V exp [0. 10(A - 12) + 3(S mod 2) + 41} , (7) 

so that e.g. the calculation of the ground-states energies of the magnesium isotopes 
24,25,26 Mg would require 0.002, 0.064, 0.131 PFlop-yr and one would need 1.9, 73.8, 145.5 
TB to store the configurations for the calculation of matrix elements for these nuclei. In 
many cases, neutron scattering off nuclei can also be calculated making use of Liischer's 
formula that relates the continuum phase shift to the finite energy shift measured on 
the lattice. The CPU and storage requirements for neutron scattering of the isotopes 
24 ' 26 Mg would be 0.064,5.43 PFlop-yr and 73.8,5473 TB, respectively. All this work 
remains to be done but bears a lot of promise. The possibilities that would be offered 
in a future exascale era are simply breath-taking. 



6 Some final words 

The chiral effective Lagrangian of QCD offers a tool to systematically investigate the 
structure of the nucleon and the nuclear interactions with high accuracy - it appears 
that we are now on the right track to achieve what could only be dreamed about a few 
decades ago - see the discussion in Sec. 1. Nuclear physics will no longer be based on 
model-building and phenomenological approaches, although still quite a bit of work is 
ahead of us. But this work - even if hard and time-consuming - will be very rewarding 
at the end. Furthermore, strangeness nuclear physics can also be addressed within this 

* 2 I am grateful to Dean Lee for providing me with these estimates. 
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framework but here a sound data basis has yet to be established. Eventually, lattice 
QCD might also contribute significantly - but there are still quite a few obstacles to be 
overcome, see e.g [47]. It is fair to say that nuclear physics as a field, and I personally, 
owe a lot to Gerry - thank you. 
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